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HIGHLIGHTS 


►  Non-isothermal  multi-D  DMFC  with  micro-porous  layers  is  developed. 

►  This  paper  presents  four  types  of  DMFC  with  anode/cathode  MPLs. 

►  The  model  explains  interfacial  saturation  jump,  water/methanol  crossover. 

►  The  model  explains  the  non-isothermal  effect  on  heat  and  mass  transfer  in  DMFCs. 

►  Large-scale  simulation  of  28  cm2  cell  is  conducted. 


ARTICLE 


N  F  O 


A  B  S  T  R 


C  T 


Article  history: 

Received  7  October  2012 
Received  in  revised  form 
9  December  2012 
Accepted  12  December  2012 
Available  online  27  December  2012 


Keywords 

Direct  methanol  fuel  cell 


Micro-porous  layer 


A  two-phase  three-dimensional  direct  methanol  fuel  cell  (DMFC)  model  with  the  capability  of  saturation 
jump  is  developed  in  order  to  investigate  the  effect  of  micro-porous  layers  (MPL)  on  the  water/methanol 
crossover  in  a  DMFC.  It  is  found  that  hydrophobic  MPL  on  the  cathode  side  helps  cathode  catalyst  layer 
reserve  water  whereas  hydrophobic  MPL  on  the  anode  side  blocks  liquid  flow  to  the  cathode  side,  which 
effectively  reduces  water/methanol  crossover.  As  the  water  transport  in  a  DMFC  occurs  in  two-phase, 
energy  equation  considering  the  heat-pipe  effect  with  latent  heat  is  solved  together.  Calculation  result 
revealed  thermal  conductivity  of  gas  diffusion  layer  has  a  strong  influence  on  both  temperature  and 
water  transport.  Finally,  a  large-scale  simulation  of  28  cm2  cell  is  conducted  to  examine  current  density 
and  fuel  concentration  distribution. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Inside  a  direct  methanol  fuel  cell  (DMFC)  several  transport 
phenomena  related  to  electrochemical  reactions  occur  simulta¬ 
neously.  Methanol  oxidation  reaction  (MOR)  occurs  on  the  anode 
side  whereas  oxygen  reduction  reaction  (ORR)  occurs  on  the 
cathode  side.  Water  participates  in  MOR  as  a  reactant  and  it  is 
produced  as  a  byproduct  from  ORR  in  the  cathode.  Therefore, 
theoretically,  it  is  possible  to  operate  a  DMFC  without  external 
supply  of  water  [1],  Such  situation  is  called  as  ‘water  neutral  con¬ 
dition’.  This  is  very  desirable  because  user  should  carry  only  pure 
methanol  without  water,  which  significantly  reduce  the  volume 
and  weight  of  the  DMFC  system.  In  order  to  construct  such  an  ideal 
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DMFC,  it  is  important  to  know  the  fundamentals  of  transport 
phenomena  which  occur  three-dimensionally  in  a  DMFC.  Many 
experimental  studies  about  these  transport  phenomena  in  fuel  cells 
have  been  conducted.  However,  there  is  a  limitation  of  such  ex¬ 
periments  because  it  is  difficult  to  visualize  and  measure  the  spe¬ 
cies  transport  inside  porous  media  of  a  fuel  cell,  which  three- 
dimensionally  occurs  in  very  small  scale  (thickness  of  typical 
porous  media  is  around  0.3  mm).  Multi-dimensional  fuel  cell  model 
incorporated  with  computational  fluid  dynamics  (CFD)  is  a  power¬ 
ful  tool  for  the  understanding  the  fundamentals  of  transport  phe¬ 
nomena  in  a  fuel  cell. 

Multi-dimensional  DMFC  model  framework  was  established 
by  Wang  and  Wang  [2],  and  Liu  and  Wang  [3],  Wang’s  model 
was  two-dimensional  and  isothermal.  Liu's  model  extended 
two-dimensional  model  to  three-dimensional  but  it  was  still  iso¬ 
thermal  model.  Capability  of  non-isothermal  modeling  is  important 
for  predicting  non-uniform  behavior  such  as  evaporation/con¬ 
densation  of  water  and  latent  heat  effect  on  heat  removal  in  a  fuel 
cell.  In  addition,  both  of  previous  models  didn’t  consider  the  effect 
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of  micro-porous  layer  which  is  proved  to  be  very  effective  for 
minimizing  water/methanol  crossover  by  generating  saturation 
jumps  at  the  interface  of  porous  media. 

In  this  paper,  firstly  we  introduce  the  two-phase  non-iso- 
thermal  multi-dimensional  DMFC  model,  which  is  developed  to 
identify  and  simulate  the  transport  phenomena  and  to  predict  the 
cell  performance  and  the  fuel  efficiency.  Secondly,  with  the 
developed  model,  the  effect  of  micro-porous  layer  on  water/ 
methanol  crossover  is  simulated  and  visualized  by  applying  satu¬ 
ration  jump  sub-model.  Thirdly,  non-isothermal  behavior  such  as 
cathode  side  heat  removal  by  the  heat-pipe  effect  is  simulated  and 
discussed.  Finally,  a  full-scale  simulation  of  28  cm2  DMFC  cell  is 
introduced. 


2.  Physical  model 

The  present  multi-dimensional  model  is  extended  from  the 
work  by  Wang  and  Wang  [2],  and  Liu  and  Wang  [3],  Important 
extended  capabilities  are  saturation  jump,  non-isothermal  behav¬ 
ior  and  large-scale  simulation  with  parallel  processing.  Governing 
equations  for  the  present  multi-dimensional  DMFC  model  are 


summarized  in  Table  1,  and  constitutive  relationships  are  listed  in 
Table  2. 


2.1.  Mass  transport  in  a  typical  DMFC 

The  structure  of  a  single  channel  in  a  typical  DMFC  is  presented 
in  Fig.  1  which  is  created  for  the  present  multi-dimensional  model. 
On  the  anode  side,  liquid  fuel  mixture  is  introduced  from  the  inlet 
of  the  anode  channel  by  a  liquid  pump.  Introduced  liquid  fuel  is 
transported  to  the  anode  GDL  surface  and  spread  throughout  the 
anode  porous  media,  then  consumed  in  the  anode  catalyst  layer  by 
MOR.  On  the  cathode  side,  oxidizer  (typically  air)  is  introduced 
from  the  inlet  of  the  cathode  channel  by  a  blower.  Introduced 
oxidizer  is  diffused  throughout  the  cathode  porous  media  and 
consumed  in  the  cathode  catalyst  layer  by  ORR. 

2.2.  Description  of  mass  transport  on  the  anode  GDL  surface 

Several  studies  [2,4,5]  explain  mass  transfer  at  the  anode 
channel/GDL  interface  by  convective  mass  transfer  which  can  be 
obtained  by  heat  transfer  analogy  as  follows: 


m 1  =  hmJ(qA  -  ciB)  =  Sh'°'  (c,A  -  C/g)  (Z.H.  Wang,  2003) 

<  =  hm  i(pi  -  ft>)  =  ShD^mJ  (Pi  -  ft>)  (J.H.  Nam,  2003) 

rh^o  =  ($2°  -  Cg2°)  =  Shf>R —  (c^2°  -  Cg2°)  (U.  Pasaogullari,  2004) 


/(Re) 

2.693  (Fully  developed  laminar  flow) 


(1) 


Source  terms 

Mass 

V-(pu)  =  S™ 

5™  =MMe0HsME0H+MH20s„20+^ 

Sg,  =  MH2°SH2°  +  M°2S°2  +MC°2J™r 

Momentum 

V./plTt?)  =  -Vp  +  V-^  +  S” 

S“  =  -|t t  (in  porous  media  domain) 

ater  transport 

V-(pi?yHj0)  =  V-[pDcopi„vyHi0]+SH2° 

Methanol  transport 

v.(^Me°H)  =  v.[/)DMeOHv.yMeOH]+sMeOH 

S“f°»  =  -jpMMeOH 

Oxygen  transport 

v-(PTtY°D  =  v-Ip^vy^i+v.^Jl  +s°* 

^  -  ->“■  -JTM°2 

Proton 

o  =  V-(K','ffVM+St 

<00  =i 

<cd  =  -Jc+jxover 

Electron 

0  =  V.(ks,eJV«+S? 

St,acl  =  ~} 

<cc,  -k-Uer 

Energy 

v-(pCpi?r)  =  v-(k^vr)  +  sT 

where  {  **  =  ^  + 

**=j(v  +  T%£)  + ^ 

sr 

^b,r,0=j(v  +  T^)  +^+j*™r(v  +  T% f) 

SU  = 

S.  Jung  /  Journal  of  Power  Sources  231  (2013)  60—81 


Parameters 

Expression 

Reference 

Relative  permeability 

K,t  -  s";  krjr  =  (1  -S) 

[15] 

Methanol  diffusivity  (vapor) 

dM0OH,HjO  =  DMeOH,C02  = 

,  „„  .r  T  \18231.013  X  105Pa  ,  2(, 

=  196  X  10  (328.15)  p  [m  /Sl 

[16] 

Water  diffusivity  (vapor) 

OjWO,  =  2Q1  x  10-5^ 

T  \18231.013  x  105Pa  r  2  ,  , 

307)  P 

[17] 

Methanol  diffusivity  (liquid) 

DMeOH,H20  =  1  4  x  IQ-9  | 

^647.3- 298.1 5\  6  ,  7  1 

K  647.3 -T  )  I™  /5! 

[18] 

Oxygen  diffusivity  (gas) 

D°2  =  3.57  x  lO-5^ 

p™,  [m2/s] 

[17] 

Water  content  in  Nafion-membrane 

[19,20] 

Water  diffusivity  in  Nafion-membrane 
Methanol  diffusivity  in  Nafion-membrane 
Water  EOD  coefficient  in  Nafion-membrane 


Water  reference  EOD  coefficient  in  Nafion 
Methanol  drag  coefficient  in  Nafion-membrane 


22  (s  >  0.3) 

{  14  +  8S/0.3  (s  <  0.3) 

l  0.043  +  17.81RH  -  39.85 RH2  +  36.0 RH3  (vapor) 


=  4.80  x  lO-i'exp^lS^ -1)]  [m2/s] 

D“f°H  =  1.5  x  l0-'°exp [2416Q3 -  m2/s] 

n«2°  =  (  [^]  (agj  -  l)  + 1  <for  ^  14> 

[l  (for  7  <  14) 

} 

nd\°f  =  1  6767  +  0.0155(r  -  273)  +  8.9074  x  10-5(T  -  273)2 

M  OH  H  0CMe°H 

km  =  0.1  (S/cm)  @  60° C  (fully  hydrated) 

Dgljy  =  Dj?2°(l  -  s)V 

1.263s3  +  1.669s2  -0.966s +  0.56  (when  flc  <  90' ) 
\  1.263s3  -  2.120s2  +  1.417s  (when  8C  >  90‘) 


m  =  \ 


However,  the  basic  assumption  of  this  description  is  ‘permeable 
surface’.  There  is  no  information  about  hydrophobicity  or  perme¬ 
ability  in  this  description.  In  other  words,  this  expression  works 
well  when  GDL  surface  has  no  hydrophobicity.  However,  many 
experimental  results  show  good  agreement  with  this  description 
because  actual  GDL  surface  is  not  ideally  hydrophobic  or 


hydrophilic  since  it  is  a  mixture  of  carbon  (hydrophilic)  and  PTFE 
(hydrophobic).  If  GDL  surface  is  ideally  hydrophobic  (liquid 
impermeable),  this  description  fails. 

Another  way  to  explain  the  mass  transport  is  applying  hydraulic 
pressure  acting  on  the  GDL  surface.  Displacing  gas-phase  by  liquid- 
phase  in  hydrophobic  porous  media  is  called  as  drainage.  There 
exist  three  methods  to  cause  drainage:  1 )  gravitation,  2)  centrifugal 
force,  3)  externally  applied  liquid  pressure.  In  a  fuel  cell,  first  and 
second  methods  are  not  applicable.  In  order  to  fill  ideally  hydro- 
phobic  porous  media  with  liquid,  we  should  apply  high  liquid 
pressure  externally  which  is  called  ‘break-through  pressure’,  which 
is  defined  as  follows: 


Breakthrough  pressure  of  GDL  can  be  estimated, 


s^c,gdl  _  4  X  0.06  X  cos!10° 


(3) 


This  is  pretty  high  value  compared  to  the  pressure  applied  in  the 
anode  channel  of  a  conventional  DMFC  (few  Pa).  In  addition,  there 
exists  CO2  gas  pressure  in  the  anode.  If  external  liquid  pressure  is 
applied  to  permeate  the  GDL  surface,  it  should  be  at  least, 


Pl>Pb  +  Pg  (4) 

But,  actual  mass  transfer  at  the  channel/GDL  interface  occurs 
without  such  a  high  external  pressure.  Therefore,  surface  liquid 
transport  of  a  conventional  DMFC  can  be  explained  by  convective 
mass  transfer  description.  If  GDL  surface  is  ideally  hydrophobic 
(convective  mass  transfer  fails)  and  channel  is  not  pressurized 
enough  (less  than  the  breakthrough  pressure  of  the  GDL  surface), 
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then  liquid  transport  is  theoretically  impossible  and  anode  will  not 
be  wet. 

On  the  cathode,  as  water  is  continuously  produced  by  ORR, 
liquid  water  should  be  removed  from  the  reaction  site.  Otherwise, 
liquid  water  will  block  the  reaction  site  and  ORR  will  be  stopped. 
Liquid  water  removal  from  the  cathode  catalyst  layer  to  the  surface 
of  the  cathode  channel  can  be  easily  explained  by  hydraulic  pres¬ 
sure  built  by  water  production  by  ORR.  As  it  is  difficult  for  liquid 
water  to  permeate  proton  exchange  membrane  which  has  very  low 
hydraulic  permeability  (or  very  high  breakthrough  pressure),  liquid 
water  produced  in  the  cathode  catalyst  layer  must  move  toward 
cathode  GDL  which  has  relatively  lower  breakthrough  pressure.  In 
addition,  as  the  applied  gas  (air)  pressure  in  the  cathode  is  almost 
constant  (~  ambient  pressure),  liquid  water  can  easily  move 
through  the  hydrophobic  GDL  by  drainage  explained  previously. 

2.3.  Description  of  mass  transport  in  the  anode  porous  media 


It  is  difficult  to  explain  liquid  transport  in  the  anode  of  a  con¬ 
ventional  DMFC  with  drainage  theory  because  of  following  reasons. 
First,  there  is  no  or  very  low  externally  applied  liquid  pressure  in  the 
anode  channel.  Second,  during  liquid  drainage  in  hydrophobic 
porous  media,  liquid-phase  pressure  must  be  greater  than  gas-phase 
pressure.  However,  in  the  anode  of  a  conventional  DMFC,  gas-phase 
pressure  (CO2  pressure)  must  be  greater  than  the  liquid  pressure  in 
the  channel.  Otherwise,  CO2  gas  cannot  be  removed  out  of  the  cell. 
Now  we  introduce  phase-diffusion  theory  (M2  model)  developed  by 
Wang  and  Cheng  [6]  in  order  to  explain  liquid  transport  in  the  anode 
of  a  DMFC  here.  M2  model  is  appropriate  for  explaining  multi-phase 
transport  occurring  in  porous  media,  it  is  widely  used  for  fuel  cell 
modeling  framework  [2,3,5,10,14,21  ].  M2  model  considers  the  multi¬ 
phase  system  as  a  chemical  mixture.  Based  upon  this,  M2  model 
describes  the  multi-phase  flow  in  terms  of  mass-averaged  mixture 
flux  velocity  and  phase-diffusive  flux  as  follows: 

Pl^l  =  7i*kpt 

PgUg  =  7g  +  hPU  (5) 

,  A,AgJC  ,  -t* 

where  j  (  =  s  Vpc  and  j  ,  +  j  g  =  0 

Phase-diffusive  flux  (  j  ;)  is  driven  by  the  gradient  of  capillary 
pressure  (pc)  defined  as  follows: 

Pc  =  Pg~  Pi  =  <rcos  Oc(£)J(s)  (6) 

Leverette’s  ./-function  (see  Table  2)  empirically  expresses  that 
capillary  pressure  is  a  function  of  liquid  saturation  defined  as 
follows: 


-=  1-- 


PgY™-PgY^ 


(7) 


Vpore  p,Y^°  -  p,yH20  +  pgyH20  _  pgY£2° 

where  two-phase  mass  fraction  relationships  are  defined, 
pYHj0  =  p,Y"2°s  +  pgYg'2°(J  -  s)  (8) 

A;  and  /.g  express  relative  mobility  of  liquid  and  gas,  respectively. 


A,  = 


n 


_  K0i 

Kj/Vt  +  hg/Vg 

(K,  ,  fer  A  1 


and  Ap-  —  1  —  Ai 


Relative  permeabilities  of  each  phase  are  expressed  by  semi- 
empirical  Bruggeman  relationship  with  liquid  saturation  as 
follows: 

krj  =  sn  and  krjr  m  (1  -  s)n  (10) 

Phase-diffusion  says  that  each  phase  moves  in  the  opposite 
direction  with  the  same  mass  flux.  In  order  to  drive  this  phase- 
diffusion,  one  of  two-phases  must  become  driving  phase  and  the 
other  phase  becomes  induced  phase.  In  the  anode,  CO2  gas  is 
produced  according  to  the  reaction  rate  (electrical  current). 
Therefore,  gas  pressure  at  A  is  higher  than  that  of  B  and  CO2  gas 
flows  from  A  to  B  in  Fig.  2.  As  liquid-phase  moves  in  the  opposite 
direction  of  the  gas  flow  by  the  definition  of  phase-diffusion  of  M2 
model,  liquid  phase  induced  from  the  anode  channel  moves  from 
B  to  A. 

CH3OH(Z)  +  H20(J)  ^C02(g)  +  6H+  +  6e 
m,  =  32  +  18  =  50g/s(B->A) 

rhg  =  -48g/s(B<-A)  (  ; 

met  =  50-48  =  2g/s(B—rA) 

When  we  consider  above  mass  flux  balance  in  the  anode  by 
assuming  the  net  water  transfer  coefficient  (a)  is  zero,  we  can  see 
the  amount  of  liquid  flow  and  gas  flow  are  similar  and  in  the 
opposite  direction.  This  supports  the  idea  that  phase-diffusion  is 
the  dominant  mass  transport  mechanism  in  the  porous  media. 

2.4.  Two-phase  water  transport  equation 

As  typical  DMFCs  operate  under  the  water  boiling  point  (100  °C), 
water  exists  and  moves  in  the  form  of  two-phase  (liquid  water  and 
water  vapor)  in  a  DMFC.  In  order  to  explain  two-phase  mass 
transfer,  we  adopt  M2  model  [6]  as  the  base  of  the  present  multi¬ 
dimensional  model  as  follows: 

v-  (PuYfc)  =  v-  [P,Diff^+pgDkgeff^]  -v-  [(yf-ygfc)T,]  +sfc 

(12) 

Left  term  represents  convectional  term  which  is  significant  in 
the  flow  channel  region.  First  term  on  the  right  hand  side  repre¬ 
sents  the  molecular  diffusive  term  and  second  term  is  phase- 


ABP 

Anode  channel 

where 


Fig.  2.  Schematic  of  phase-diffusion  in  the  anode. 


S.  Jung  /  Journal  of  Power  Sources  231  (2013)  60—81 


diffusive  term  or  capillary  flux  which  is  significant  in  the  porous 
region.  Capillary  flux  of  liquid  phase  can  be  explained  with  Darcy’s 
law  as  follows: 


(13) 


By  plugging  Eq.  (6)  into  Eq.  (13).  capillary  flux  of  liquid  phase  can 
be  expressed  as  a  function  of  liquid  saturation  as  follows: 


2.6.  MOR  and  ORR 

Butler-Volmer  equation  is  simplified  as  Tafel  equation  in  this 
study.  It  is  known  that  methanol  oxidation  reaction  (MOR)  is  Oth 
order  when  the  methanol  concentration  is  high  enough  and  1st 
order  when  the  methanol  concentration  approaches  the  depletion 
limit.  Following  equation  is  used  to  simulate  a  smooth  transition 
from  Oth  order  kinetics  to  1st  order  kinetics  in  this  study  [7-9], 


7,.^vpt  =  iMv[«os(«(l),/27(s)] 

_  AgAjffCOS (dc)(Ke)'/2  d 


(14) 


cMe°H+ifcexp^^ 


[A/m3] 


where  r/a  =  &s  -  &e  -  lf¥J 


(20) 


By  plugging  Eq.  (7)  into  Eq.  (14),  capillary  diffusive  flux  can  be 
reformed  as  follows: 


On  the  cathode  side,  oxygen  reduction  reaction  (ORR)  is 
assumed  to  be  proportional  to  oxygen  concentration. 


yt  ,=  Ag7;ffCOS(flc)(Kc)1/2  dj  Pi(PgY"2°  PgYg2°) _ vyH2( 

ds  (p,Y»>°  -  P(yH2o  +  PgyH2o  _  PgYy°]2 


We  can  assume  the  water  concentration  in  the  gas  phase  is 
always  saturated  and  the  gradient  of  mass  fraction  in  the  liquid 
phase  is  almost  zero  due  to  dilute  solution.  Therefore,  the  water 
transport  equation  can  be  derived  from  Eq.  (12)  with  Eq.  (15)  as 
follows: 

V  •  (PTtYH*°)  =  V  •  [PDcapillV  Y*°]  +  SH2°  (16) 

This  equation  explains  the  water  is  transported  by  capillary 
diffusion  in  the  porous  media  with  the  capillary  diffusivity  being 


M  (2,) 

where  pc  =  <ps  -  <I>e  -  U°  [V] 

2.7.  Methanol  crossover  rate  (MCO) 

Some  amount  of  methanol  that  couldn’t  react  in  the  anode  catalyst 
layer  may  move  to  the  cathode  side  by  permeating  the  membrane. 
This  crossover  methanol  can  participate  either  in  the  electrochemical 


DcaPill~l(  Y”*°-YW) 


AgV|cos(flc)|(fre)1/2  dj  Pi(PgYl'1°  ~  %yg2°) _ 

%yr  -  P,Y ™  +  Pg^o  -  PgYg  2°)  2 


(17) 


2.5.  Two-phase  methanol  transport  equation 

By  assuming  methanol  vapor-liquid  equilibrium  in  the  anode 
(Henry’s  law),  two-phase  mass  fraction  relationship  (Eq.  (8))  of 
methanol  can  be  reformed  as  follows: 


By  plugging  Eq.  (18)  into  Eq.  (12),  the  methanol  transport 
equation  becomes 


reaction  which  produces  crossover  current  iX0Ver  or  in  the  pure 
chemical  reaction,  i.e.  combustion.  Crossover  current  has  negative 
effect  on  the  cell  performance  since  it  causes  mixed  potential  which 
degrades  the  cell  performance.  In  addition,  crossover  current  con¬ 
sumes  additional  fuel  and  oxidizer,  which  reduces  the  fuel  efficiency 
of  DMFC  also.  Crossover  current  can  be  described  as  shown  in  Eq.  (22) 
[3  ].  It  is  considered  that  methanol  electro-osmosis  drag  and  methanol 
diffusion  contribute  to  the  total  methanol  crossover. 


Wr  =  6F<|OH 

j  :  cMeOH  (22) 

Where  N^eOH  =  W  =  nMeOH>  +  DMeOHC0d 

6F  “  F  omem 

Methanol  crossover  rate  (MCO)  describes  the  fuel  efficiency  of 
a  DMFC  as  follows  [12): 


SJung/J 


2.8.  Cell  voltage 

Generally,  cell  voltage  is  described  as  following  equation  con¬ 
sidering  overpotentials. 

Vcell  =  (u°c  -  17c)  -  (u°a  -  17  a)  ~  IRmem  ~  I 'Rather  (24) 

However,  as  the  present  DMFC  mode  simultaneously  solves  all 
the  governing  equations  listed  in  Table  1,  cell  voltage  can  be  sensed 
by  the  electron-phase  potential  value  difference  between  anode 
bipolar  plate  and  cathode  bipolar  plate  as  follows: 

Vcell  =  ®s,CBP  -  &SABP  (25) 

This  cell  voltage  includes  all  of  overpotentials  such  as  ohmic 
loss,  concentration  loss  and  kinetics  loss. 

2.9.  Net  water  transfer  coefficient  (a) 

Water  is  a  necessary  reactant  in  MOR.  Water  transport  through 
membrane  can  be  expressed  using  the  net  water  transport  coeffi¬ 
cient  defined  in  Eq.  (26)  which  consists  of  electro-osmosis  drag, 
diffusion  and  hydraulic  permeation  [10],  By  examining  the  water 
transport  coefficient,  we  can  find  out  whether  water  is  being  lost  or 
gained. 


2.10.  Multi-dimensional  saturation  jump  model 

Qi  and  Kaufman  [11]  experimentally  showed  the  effect  of 
micro-porous  layer  (MPL)  which  is  inserted  between  GDL  and 
catalyst  layer  on  the  water  transport  and  distribution  in  fuel  cells. 
As  this  MPL  is  made  of  finer  and  more  hydrophobic  materials  than 
other  porous  layers,  a  discontinuity  in  the  liquid  saturation  across 
the  interface  of  two  layers  occurs  while  the  liquid  pressure  at  the 
interface  is  continuous.  Anode  MPL  has  the  effect  of  reducing 
saturation  level  in  the  anode  catalyst  layer  by  blocking  water  flow 
from  the  anode  channel  whereas  cathode  MPL  has  the  effect  of 
increasing  saturation  level  in  the  cathode  catalyst  layer  by 
reserving  produced  water  by  ORR.  Therefore,  the  net  water  flux 
from  the  anode  to  cathode  is  greatly  reduced  or  even  it  negatively 
flows  from  the  cathode  to  anode,  which  is  helpful  for  water 
management.  Nam  and  Kaviany  [4]  theoretically  explained  the 
saturation  discontinuity  effect  with  one-dimensional  mathemati¬ 
cal  model,  so  called  saturation  jump  model.  Passaogullari  and 
Wang  [5]  further  investigated  the  effect  of  saturation  jump  with 
MPL  on  the  water  transport  in  the  porous  media  of  a  fuel  cell. 
Shaffer  and  Wang  [12]  conducted  a  theoretical  discussion  about 
the  effect  of  hydrophobic  MPL  on  the  cell  performance  and  the 
fuel  efficiency  of  DMFCs  by  simulating  one-dimensional  model.  In 
order  to  minimize  water/methanol  crossover  in  a  DMFC,  most  of 
DMFCs  recently  developed  have  two  MPLs  with  thin  proton 


N±°J-  = 


(26) 
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exchange  membrane,  so  called  ‘low-a  MEA’.  Experimental  results 
by  Liu  and  Wang  [13]  showed  that  the  cell  performance  and  fuel 
efficiency  can  be  greatly  improved  with  properly  selected  MPLs  in 
the  anode  and  cathode.  As  Nam  and  Kaviany  [4]  and  Shaffer  and 
Wang  [12]’s  one-dimensional  model  directly  solved  flux-based 
liquid  saturation  equation  with  above  relationship,  their  model 
easily  shows  a  saturation  profile  in  through-plane  direction. 
However,  as  the  present  multi-dimensional  model  solves  mass 
fraction  equation  instead  of  saturation  equation,  saturation  dis¬ 
tribution  is  acquired  from  mass  fraction  distribution  in  Eq.  (7). 
Furthermore,  several  parameters  and  constitutive  relationships  of 
the  multi-dimensional  model  depend  on  saturation  distribution, 
which  means  saturation  and  mass  fractions  are  strongly  coupled 
in  the  model.  We  introduce  how  to  implement  saturation  jump  in 
the  present  multi-dimensional  DMFC  model.  Fig.  3(a)  shows  the 
interface  of  two  different  porous  layers  (GDL  and  MPL).  As  the 
capillary  pressure  should  be  continuous  at  this  interface  as 
follows: 


ffcos^i  (^);(s1Jnt)  =  <rCOS0c,2(|)j(S2,int)  (27) 

When  the  saturation  of  one  layer  is  given,  the  saturation  of  the 
other  layer  can  be  found  by  solving  following  equation.  Calculation 
result  is  shown  in  Fig.  3(b). 


J(«  bin 


1J(S2,  ini 


01  =  ffC0S(6ll)(|^) 

h  =  <rc°s(02  )(|0  ' 


(28) 


However,  the  multi-dimensional  model  calculates  liquid  satu¬ 
ration  in  each  mesh  (si  and  S2)  and  interfacial  liquid  saturation  ( Sj„t ) 
is  not  given.  Therefore,  interfacial  liquid  saturation  should  be 
determined  by  Si  and  S2  by  water  flux  continuity  (liquid  phase  and 
vapor  phase)  at  the  mesh  interface  as  follows: 


Nint  =  N1)in  t  =  N2,in 

n  sUnt  — 51  n 


Fig.  4.  Schematic  of  heat-pipe  effect  in  a  PEMFC  [14], 


Then  from  Eq.  (29),  we  have, 
S2,int  =  ■A '■$  ijnt  +  B 


(30) 


where 


A 

B 


Ds,  1  Ax2 
~DS  2 


(d»et 


A*2 

Ds,  2 


fg,i  ~  °g,  1  S  ‘"AXl  4 


Now  we  combine  Eq.  (28)  and  Eq.  (30)  to  get  interfacial  liquid 
saturation  and  interfacial  liquid  flux.  Interfacial  liquid  saturation 
does  not  explicitly  show  up  in  the  multi-dimensional  model. 
Instead,  we  add  or  subtract  interfacial  liquid  flux  found  above  in  the 
mesh  1  and  2  as  follows: 


S,  =  -V- (pDCQpi(/VyH2°)  +M^°V-Nint 
S2  =  -V-(pDCQpi„vy^°)  -M^°VWint 


Finally,  above  water  source  term  is  put  to  the  water  equation  of 
Eq.  (16). 


2.11.  Two-phase  heat  transfer  model 

As  liquid  water  and  water  vapor  coexist  in  a  DMFC,  water  con¬ 
centration  is  strongly  affected  by  temperature  whereas  tempera¬ 
ture  is  affected  by  the  latent  heat  due  to  the  water  phase  change 
(condensation  and  evaporation).  Therefore,  species  equation  and 
energy  equation  are  strongly  coupled  in  the  present  multi¬ 
dimensional  model.  Standard  energy  equation  for  the  present 
DMFC  model  is  derived  from  M2  model  [6], 


Fig.  5.  Three-dimensional  mesh  of  the  single  channel  DMFC  model  for  the  present  study. 


S.  Jung  /  Journal  of  Power  Sources  231  (2013 )  60—81 


Table  3 

Cell  geometry  and  material  properties  of  the  baseline  DMFC  model. 


Table  4 

Simulation  parameters  and  operating  conditions  of  the  baseline  DMFC  model. 


Description 


Anode  GDL  thickness 
Anode  MPL  thickness 
Anode  CL  thickness 
Membrane  thickness 
Cathode  CL  thickness 
Cathode  MPL  thickness 
Cathode  GDL  thickness 
GDL  porosity 
MPL  porosity 
CL  porosity 
FD  porosity 
GDL  permeability 
MPL  permeability 
CL  permeability 
FD  permeability 
Membrane  permeability 
GDL  contact  angel 
MPL  contact  angel 

Channel  width 
Channel  height 
Land  width 
Bipolar  plate  thickness 


Kgdl 

«mpl 

Kd 

Kfd 

Kme„ 


Lcell 

Warn, 

Wland 

HbP 


270  pm 
30  pm 

45  pm 
25  pm 
30  pm 
210  pm 
0.6 
0.4 
0.4 
03 

2.0  x  1(T12  m2 
5.0  x  1(T14  m2 
1.0  x  1(T13  m2 
2.84  x  10-13  m2 


114” 

96” 

75  mm 


0.5  mm 


v-  (pcpiiT)  =  v-(/<vr)  +  sT 

where  ST  =  STsensible  +  Sjatent 


(32) 


There  are  three  sensible  heat  source  terms  in  the  anode:  1) 
irreversible  heat  of  the  electrochemical  reaction,  2)  reversible 
entropic  heat,  3)  ohmic  heating,  respectively  (see  Eq.  (33)).  In  the 
membrane,  only  ohmic  heating  exists  since  electrochemical  reac¬ 
tion  does  not  occur  there. 


Sensible, act  =J'(??  +  T  dT  ) 
c T 

spmihlp  mpi 


T6S)+^ 


tv*  -j(’ 


(33) 


The  influence  of  methanol  crossover  on  energy  balance  is  added 
also,  which  generates  additional  heat  in  the  cathode  catalyst  layer 
where  crossovered  methanol  reacts.  Latent  heat  effect  is  included 
in  the  present  model  as  following  equation. 

^latent  =  Vfg  = 

=  V-  ^VC0S(gc)(ye)1/2gVs  _  hy£vp) 

(where  riifg  is  the  phase  -  changing  ratej  (34) 


Reference  current  density 

Anode  reference  exchange  current  density 

Cathode  reference  exchange  current  density 

Anode  catalyst  ionomer 

Cathode  catalyst  ionomer 

Anode  catalyst  loading 

Cathode  catalyst  loading 

Surface  tension 

Inlet  methanol  concentration 

Pure  liquid  methanol  concentration 

Pure  liquid  water  concentration 

Total  anode  transfer  coefficients 

Total  cathode  transfer  coefficients 

Faraday  constant 

Universal  gas  constant 

Equivalent  weight  of  dry  membrane 

Dry  membrane  density 

Cell  temperature 

Anode  stoichiometry 

Cathode  stoichiometry 

Operating  current  density 


CMeOH.pure 

CH20.pure 


150  mA  cm~2 
45.51  A  m2 
0.12  Am“2 
25  vol.% 

25  vol.% 

4  mg  cm-2 
2  mg  cm  2 
0.0625  N  m-1 
4000  mol  nr3 
24,719  mol  m3 
55,556  mol  m-3 
0.239 
0.875 

96,487  C  mor1 
8.314  J  mol-1  K 
1.1  kg  mor1 
1980  kg  m-3 
313  K 

2.0 

150  mA  cm“2 


This  source  term  arises  from  latent  heat  release  or  absorption 
due  to  water  condensation/evaporation  in  two-phase  flow  region. 
Evaporation  occurs  in  the  near  surface  area  of  GDL  where  air  flow 
carries  moist  away.  Condensation  occurs  in  the  cold  region  (under 
the  land  region)  and  cathode  catalyst  layer  region  where  water  is 
being  produced  by  ORR.  Because  the  present  model  assumes  water 
is  produced  first  in  the  vapor  form,  excessive  water  which  is  greater 
than  the  saturation  value  must  be  condensed  as  liquid  form. 

2.12.  Non-isothermal  behavior:  the  heat-pipe  effect 

In  a  typical  DMFC,  heat  transfer  in  the  anode  is  mostly  carried  out 
by  liquid  flow  in  the  channel  which  has  large  convective  heat  transfer 
coefficient  (liquid  cooling).  However,  in  the  cathode,  heat  conduction 
and  heat-pipe  effect  are  important  since  there  is  no  liquid  flow. 
Basu  et  al.  [14]  investigated  the  phase  change  and  heat-pipe  effect 
in  a  hydrogen  PEMFC  (proton  exchange  membrane  fuel  cell).  They 
showed  the  cell  operating  condition  ( relative  humidity  of  the  inlet  gas) 
and  material  property  (heat  conductivity  of  GDL)  have  a  strong  influ¬ 
ence  on  the  phase  change  and  temperature  distribution  by  quantifying 
the  phase-changing  rate  and  the  portion  of  heat  transfer  by  heat-pipe 
effect.  In  this  part,  we  reproduce  the  non-isothermal  behavior  with 
the  heat-pipe  effect  in  a  DMFC  by  adopting  Basu’s  methodology. 

Fig.  4  shows  how  heat  transfer  occurs  in  the  cathode  by  heat- 
pipe  effect.  As  the  catalyst  layer  is  hottest  due  to  reaction  and 
land  is  coldest,  major  heat  flow  direction  is  from  the  catalyst  to 
GDL/land  interface.  Water  vapor  near  the  catalyst  layer  is  trans¬ 
ported  to  GDL/land  interface  by  thermal  diffusion  and  it  condenses 
there  due  to  lower  temperature.  Condensed  water  near  GDL/land 
interface  moves  toward  the  gas  channel  and  finally  it  evaporates 
since  dry  air  flow  carries  moist  away.  By  considering  heat  balance 
near  GDL/land  interface,  following  result  can  be  derived: 


*  [(v. + '4""-0 £]  -  *  ("*£) 

!Dg  eff  =  Dg.2°(l  -  s)2  3e4  :  Effective  water  vapor  diffusivity  (35) 

keff  =  kgdi  +  hfgMH2°Dg2eff  '■  Effective  thermal  conductivity 
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Above  result  implies  that  the  heat-pipe  effect  increases  the 
effective  thermal  conductivity,  which  improves  heat  release. 

In  addition  to  the  heat  transfer  effect,  heat-pipe  effect  includes 
water  transport  effect  also.  When  the  cell  is  non-isothermal,  ther¬ 
mal  vapor  diffusion  occurs  due  to  water  vapor  concentration 


gradient  in  the  two-phase  region.  As  this  thermal  vapor  diffusion 
occurs  from  hot  to  cold  region,  liquid  saturation  in  the  hot  region 
decreases  and  it  increases  in  the  cold  region,  which  means  water 
removal  from  the  cathode  porous  media  to  the  cathode  channel  is 
enhanced.  As  water  is  transported  by  both  capillary  diffusion  and 


Fig.  7.  Methanol  concentration  in  the  anode  catalyst  layer  of  the  single  channel  cell  model  according  to  the  flow  stoichiometry:  (left)  Tceu  =  313  K,  '(a  =  1.4,  'C<:  =  2.0,  (center) 
Tcell  =  313  K,  =  2.4,  =  2.4,  (right)  Tce„  =  313  K,  =  5.0,  ?c  =  5.0. 
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Fig.  8.  Crossover  current  density  distribution  of  the  single  channel  cell  model  according  to  the  flow  stoichiometry:  (left)  Tceii  =  313  K,  =  1.4,  £c  =  2.0,  (center)  Tceii  =  313  K,  £a 
Cc  =  2.4,  (right)  T«„  =  313  K,  =  5.0,  =  5.0. 


Fig.  9.  Current  density  distribution  of  the  single  channel  cell  model  according  to  the  flow  stoichiometry:  (left)  Ten  =  313  K,  =  1.4,  ?c  =  2.0,  (center)  Tce„  =  313  K,  =  2.4, 
(right)  Tce„  =  313  K,  =  5.0,  ?c  =  5.0. 
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non-isothermal  vapor  diffusion  in  the  two-phase  region,  we  can 
expect  the  net  water  transfer  coefficient  will  increase.  Contrary  to 
the  heat  equation,  those  two  water  transport  mechanism  (capillary 
diffusion  and  thermal  vapor  diffusion)  are  not  compatible.  Hence, 
they  cannot  be  incorporated  in  a  single  diffusion  term.  Therefore, 
we  introduce  an  additional  water  source  term  to  reflect  the  thermal 
vapor  diffusion  as  follows: 

sm  =  /v'  [CCm]  (36) 

This  term  is  added  to  the  water  equation  Eq.  (16). 


2.13.  Full-scale  DMFC  model  with  parallel  channel  configuration 

Fig.  5  shows  the  mesh  used  for  the  present  multi-dimensional 
DMFC  modeling,  which  simulates  a  single  channel  (75  mm  length) 
of  a  28  cm2  DMFC  which  has  18  parallel  channels.  In  order  to  save 
computational  time,  we  assume  every  single  channel  is  identical  in 
the  parallel  cell  configuration  without  any  channel  clogging  by  CO2 
bubbles.  Therefore,  this  single  channel  can  be  considered  to  represent 
the  entire  DMFC  cell  area.  Cell  specification  and  simulation  parame¬ 
ters  are  listed  in  Tables  3  and  4,  respectively. 


2.14.  Full-scale  DMFC  model  with  serpentine  channel  configuration 

Many  DMFCs  adopt  serpentine  channel  configuration  in  order  to 
avoid  channel  clogging  problem  although  serpentine  channel 
causes  larger  pressure  drop  than  the  parallel  configuration.  On  the 
cathode  side,  removed  water  from  the  GDL  surface  may  condense 
in  the  cathode  channel,  which  may  clog  the  channel  and  hinder  air 
supply  although  it  is  not  so  much  severe  as  the  anode  side  since 
a  typical  DMFC  operates  at  low  current  compared  to  hydrogen 
PEMFC  (water  production  is  almost  1/10).  Therefore,  DMFC  with 
serpentine  channel  configuration  will  have  non-uniform  distribu¬ 
tion  of  species,  current  and  temperature  compared  to  DMFC  with 
parallel  channel  configuration. 

Fig.  6  shows  3-D  mesh  of  the  full-scale  DMFC  model  (28  cm2). 
Anode  has  single-pass  serpentine  channel  in  order  to  sweep  away 
CO2  bubbles  whereas  cathode  has  two-pass  serpentine  channel.  As 
volumetric  air  flow  rate  is  much  larger  than  liquid  fuel  flow  rate, 
pressure  drop  in  the  cathode  channel  is  much  greater  than  that  of  the 
anode  channel.  In  order  to  reduce  pressure  drop  due  to  this  large  air 
flow  rate,  2—4  serpentine  passes  are  arranged  to  grasp  both  reduced 
pressure  drop  and  better  removal  of  water  droplets.  Anode  channel 
and  cathode  channel  are  arranged  as  cross-flow  configuration  in  or¬ 
der  to  minimize  non-reactive  area  and  promote  uniform  distribution 
of  fuel  and  oxygen. 


Fig.  10.  Current  density  distribution  in  the  cf 
He  =  2.4,  (c)  rce„  =  313  K,  (a  =  5.0,  He  =  5.0. 


1  parallel  channel  configuratic 


1  the  flow ! 


liometry:  (a)  =  313  K,  'Ca  =  1.4,  =  2.0,  (b)  =  313  K,  '(a  =  2.4, 
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Table  5 

Effect  of  flow  stoichiometry  on  the  cell  performance  (single  channel  model). 

Tceii  (K)  RH  CMeoH  Zc  i  (mA  Vceii  MCO 

(*)  (M)  cm'2)  (V)  (%) 

Case  A  313.0  10.0  4.0  1.4  2.0  150  0.409  12.0 

Case  B  313.0  10.0  4.0  2.4  2.4  150  0.455  19.1 

CaseC  313.0  10.0  4.0  5.0  5.0  150  0.461  28.1 


2.15.  Implementation 

Present  mathematical  model  is  implemented  into  a  commercial 
CFD  solver  STAR-CD™  with  user  subroutine  capability.  It  is  consid¬ 
ered  the  convergence  is  achieved  when  the  residuals  reach  1CT6. 
Parallel  computing  is  used  for  the  full-scale  simulation. 

3.  Results  and  discussion 

3.1.  Parallel  channel  configuration 

3.1.1.  Effect  of  flow  stoichiometry  on  the  cell  performance  and  fuel 
efficiency 

As  the  flow  stoichiometry  (or  flow-rate)  increases,  methanol 
concentration  in  the  anode  increases  also  as  shown  in  Fig.  7.  Under 
the  channel  region  has  higher  methanol  concentration  than  under 


the  land  region  due  to  the  short  diffusion  length.  As  defined  in  Eq. 
(22),  crossover  current  density  is  strongly  affected  by  the  methanol 
concentration  near  membrane/ACL  interface.  Therefore,  crossover 
current  density  distribution  is  similar  to  that  of  methanol  concen¬ 
tration  as  shown  in  Fig.  8. 

Fig.  9  shows  the  current  density  distribution  in  a  single  channel 
model  according  to  the  fuel  stoichiometry.  When  fuel  supply  is  large 
enough,  current  density  is  uniformly  distributed  in  the  flow  stream 
direction  as  shown  Fig.  9(c)  and  the  highest  current  density  is 
observed  under  the  land  region.  In  a  DMFC,  when  methanol  supply  is 
large  enough,  MOR  is  0th  order,  which  means  MOR  is  not  so  much 
affected  by  the  methanol  concentration  in  the  anode  catalyst  layer. 
Flowever,  excessive  methanol  transported  to  the  cathode  side  (under 
the  channel  region)  degrades  the  cell  performance  by  producing 
mixed  potential,  which  results  in  lower  current  density  under  the 
channel  region  where  methanol  crossover  is  high.  When  fuel  stoi¬ 
chiometry  is  low,  MOR  is  affected  by  the  methanol  concentration  in 
the  anode  catalyst  layer  and  the  highest  current  density  is  observed 
under  the  channel  region  as  shown  in  Fig.  9(a).  If  fuel  supply  is  in¬ 
termediate,  the  highest  current  density  shows  up  under  the  land  in 
the  inlet  region  and  it  is  observed  under  the  channel  in  the  outlet 
region  as  shown  in  Fig.  9(b).  Although  large  stoichiometry  produces 
better  cell  performance,  the  fuel  efficiency  decreases  due  to 
large  MCO.  In  addition,  operating  a  cell  with  large  stoichiometry 
requires  much  pumping  power  (parasite  loss).  Contrarily,  too  low 


Fig.  11.  Liquid  saturation  distribution  at  the 


B,  C:  (a)  no  MPL,  (b)  anode  MPL  only,  (c)  cathode  MPL  only,  (d)  anode  and  cathode  MPLs. 
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Table  6 

Effect  of  MPL  on  the  cell  performance  (single  channel  model). 


Anode  MPL  Cathode  MPL  Tce„(K)  RH(%)  CMeoH(M)  ifmAcnr2)  Vce„(V)  Alpha  MCO(%) 

Case  A  NO  NO  333.0  900  15  150.0  TO  TO  0.399  061  025 

Case  B  YES  NO  333.0  90.0  2.0  150.0  3.0  3.0  0.407  -0.08  0.20 

Case  C  NO  YES  333.0  90.0  2.0  150.0  3.0  3.0  0.400  0.55  0.25 

Case  D  YES  YES  333.0  90.0  2.0  150.0  3.0  3.0  0.407  -0.20  0.20 


stoichiometry  may  lead  to  large  local  concentration  overpotential  due 
to  anode  non-uniformity  which  significantly  degrades  the  cell  per¬ 
formance  although  the  fuel  efficiency  is  good  and  the  parasite  loss  is 
minimal. 

Finally,  Fig.  10  shows  the  current  density  distribution  on  the 
entire  cell  area  that  we  can  expect  when  every  single  channel  is 
identical.  Parallel,  strip-shaped  distribution  is  observed  due  to  the 
channel  configuration.  Again,  large  flow  stoichiometry  results  in 
more  uniform  distribution  of  current  density.  Calculation  results 
about  stoichiometry  effect  are  summarized  in  Table  5. 

3.1.2.  Liquid  saturation  distribution  and  effect  of  MPL 

Fig.  11  shows  liquid  saturation  distribution  in  the  cell  cross- 
section  according  to  the  MPL  configuration.  Anode  porous 
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media  is  gradually  getting  gaseous  in  the  channel  direction  since 
produced  CO2  gas  accumulates  whereas  cathode  porous  media 
becomes  wetter  since  produced  water  accumulates.  In  the 
anode,  the  region  under  the  channel  is  wetter  than  the  region 
under  the  land.  This  is  because  CO2  removal  is  blocked  by  the 
land  and  the  region  under  the  channel  directly  faces  liquid 
flow  in  the  channel.  Contrarily,  the  region  under  the  land  on  the 
cathode  side  is  wetter  than  the  region  under  the  land  since 
liquid  water  removal  from  the  cathode  catalyst  layer  to  the 
cathode  channel  is  blocked  there.  Therefore,  it  is  expected  that 
net  water  transfer  coefficient  will  be  positive  under  the  channel 
and  negative  (or  small)  under  the  land.  When  there  are  no  MPLs, 
two  weak  saturation  jumps  are  observed  at  CL/GDL  interface  as 
shown  in  Fig.  11(a). 


Fig.  12.  Liquid  saturation  profile  under  the  channel  at  the 


B,  C:  (a)  no  MPL,  (b)  anode  MPL  only,  (c)  cathode  MPL  only,  (d)  anode  and  cathode  MPLs. 
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When  the  anode  MPL  is  inserted  between  the  anode  GDL  and 
the  anode  catalyst  layer,  we  can  see  a  strong  saturation  jump  at 
the  interface  between  the  anode  MPL  and  the  anode  catalyst 
layer,  which  effectively  blocks  liquid  flow  from  the  anode  GDL  to 
the  anode  catalyst  layer  and  lowers  the  net  water  transfer  to  the 
cathode  (see  Fig.  11(b)).  In  addition,  the  anode  MPL  plays  as 
a  methanol  transport  barrier  which  reduces  methanol  concen¬ 
tration  in  the  anode  catalyst  layer  and  lowers  methanol 
crossover. 

Fig.  11(c)  shows  the  effect  of  the  cathode  MPL  inserted  between 
the  cathode  catalyst  layer  and  the  cathode  GDL.  By  reserving  liquid 
water  in  the  cathode  catalyst  layer,  the  cathode  MPL  builds  up 
hydraulic  pressure  and  increases  water  back-diffusion  to  the  anode 
side,  which  lowers  the  net  water  transfer  coefficient  also.  Three 
saturation  jumps  are  observed  in  this  case.  Note  that  liquid 
saturation  in  the  cathode  catalyst  layer  is  especially  increased. 
Cathode  MPL  does  not  have  the  effect  reducing  methanol  crossover 
(Table  6). 

When  a  cell  has  MPL  on  the  anode  and  the  cathode  side,  four 
saturation  jumps  can  be  observed  as  shown  in  Fig.  11(d).  Perfor¬ 
mance  comparison  according  to  the  MPL  configuration  is  summa¬ 
rized  in  Table  6.  Both  MPL  has  the  effect  to  lower  water  crossover 


but  only  anode  MPL  has  additional  effect  to  lower  methanol 
crossover.  Case  D  which  has  two  MPLs  shows  the  lowest  water / 
methanol  crossover. 

Figs.  12  and  13  show  the  effect  of  MPL  in  a  different  viewpoint  as 
a  form  of  saturation  profile.  Saturation  jumps  are  more  clearly 
noticed.  When  the  cathode  MPL  is  inserted,  liquid  saturation  level 
in  the  cathode  GDL  is  lowered  due  to  the  increased  water  back- 
diffusion.  This  may  result  in  wet-to-dry  transition  near  outlet  re¬ 
gion  as  shown  in  Fig.  12(c)  and  (d). 

3.1.3.  Effect  of  the  thermal  diffusion  and  phase-change  on  heat  and 
water  management 

Fig.  14(a)  and  (b)  compares  temperature  distribution  and 
thermal  diffusion  rate  in  the  cathode  according  to  the  heat  con¬ 
ductivity  of  the  cathode  GDL.  Low  heat  conductivity  of  GDL 
(Fig.  14(a))  results  in  higher  temperature  difference  than  high  heat 
conductivity  case  (Fig.  14(b)).  This  high  temperature  difference 
drives  much  more  water  thermal  diffusion  from  the  hot  region 
(near  the  cathode  catalyst  layer)  to  the  cold  region  (GDL/CBP 
interface  region).  Note  negative  value  denotes  water  is  being  lost 
whereas  positive  value  denotes  water  is  being  supplied  there.  The 
water  transported  to  the  region  GDL/CBP  interface  is  carried  away 


Fig.  13.  Liquid  saturation  profile  under  the  land  at  the  i 


s-section  A,  B,  C:  (a)  no  MPL,  (b)  anode  MPL  only,  (c)  cathode  MPL  only,  (d)  anode  and  cathode  MPLs. 
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Fig.  14.  Thermal  diffusion  rate  in  the  cathode:  (a)  kgd,  =  1.0  W  m"1  K.  RH  =  30%,  (b)  kgdl  =  3.0  W  m”1  K,  RH  =  30%,  (c)  kgd,  =  1.0  W  nr1  K,  RH  =  70%. 


by  the  air  flowing  through  the  cathode  channel.  As  the  higher 
temperature  difference  causes  much  more  water  transport  (higher 
net  water  transfer  coefficient  value),  it  can  be  concluded  that  low 
GDL  heat  conductivity  leads  to  higher  net  water  transfer  coefficient 
(a).  Effect  of  supply  air  humidity  on  the  thermal  diffusion  rate  is 
presented  by  Fig.  14(a)  and  (c).  When  highly  humidified  air  is 
flowing  in  the  cathode  channel  (CGC),  less  thermal  diffusion  occurs 
due  to  the  reduced  water  concentration  difference.  Therefore,  less 
water  is  transported  from  CCL  to  CGC,  which  results  in  lower  net 
water  transfer  coefficient  value  (a  =  0.58)  compared  to  the  dry  air 


case  (RH  =  10%)  that  induces  more  water  to  the  cathode  side 
(a  =  1.20). 

The  phase-change  rate  is  presented  in  Fig.  15.  Positive  value 
denotes  condensation  whereas  negative  value  denotes  evaporation. 
Water  vapor  transported  from  the  cathode  catalyst  layer  to  CBP/ 
GDL  interface  by  thermal  diffusion  condenses  and  it  moves  toward 
the  cathode  channel  and  evaporates.  Phase-change  rate  is  affected 
by  both  material  property  (heat  conductivity)  and  operating 
condition  (relative  humidity).  When  heat  conductivity  is  low 
(Fig.  15(a))  more  water  condenses  at  CBP/GDL  interface  and 
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Fig.  15.  Temperature  and  phase-change  rate  in  the  cathode:  (a)  ksd,  =  1.0  W  m-'  K,  RH  =  30%,  (b)  kgdl  =  3.0  W  m-'  K,  R H  =  30%,  (c)  ksdl  =  1.0  W  m1  K,  RH  =  70%. 


evaporates  at  CGC/GDL  interface  compared  to  the  high  heat 
conductivity  case  (Fig.  15(b))  since  more  water  is  transported  by 
increased  temperature  difference  with  low  heat  conductivity.  Effect 
of  supply  air  humidity  on  the  phase-change  rate  is  presented  by 
Fig.  15(a)  and  (c).  Increased  relative  humidity  of  the  air  flowing  in 
the  channel  results  in  higher  liquid  saturation  and  lowered  satu¬ 
ration  gradient  in  the  cathode  porous  media,  which  leads  to  low¬ 
ered  phase-change  rate  (see  Eq.  (34)).  From  this,  we  can  see  heat 
and  water  transport  are  strongly  coupled  by  both  material  property 
(heat  conductivity)  and  cell  operating  condition  (temperature  and 
humidity). 


Table  7 

Calculation  summary  of  full-scale  DFMC  simulation  with  serpentine  channels 
(28  cm2). 

Tce„(K)  RH  CMeoH  Cc  i  (mA  Vcel,  MCO 

(%)  (M)  cm"2)  (V)  (%) 

Case  A  313.0  10.0  4.0  2.0  2.0  150.0  0.459  13.7 

Case  B  313.0  10.0  4.0  1.4  2.0  150.0  0.428  9.2 

CaseC  313.0  10.0  4.0  2.0  5.0  150.0  0.466  14.4 
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3.2.  Serpentine  channel  configuration  (28  cm2  cell) 

In  order  to  see  the  effect  of  flow  stoichiometry,  three  cases  are 
simulated.  Simulation  conditions  are  listed  in  Table  7. 


3.2.1.  Methanol  concentration  and  crossover  current  distribution 
Fig.  16  shows  methanol  concentration  distribution  in  the 
anode  catalyst  layer.  Case  B  (?0  =  1.4)  shows  severe  non- 
uniform  distribution  of  methanol  concentration,  which 


a 


Fig.  16.  Methanol  concentration  (mol  nrr3)  in  the  anode  catalyst  layer  of  the  full-scale  cell  (28  cm2)  when  Tceu  =  313  K ,RH=  10%,  CMeoH  =  4  M,  and  1  =  150  mA  cm2:  (a)  =  2.0  and 

=  2.0  (case  A),  (b)  £0  =  1.4  and  ?c  =  2.0  (case  B),  (c)  =  2.0  and  ?c  =  5.0  (case  C). 
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implies  that  the  outlet  region  may  suffer  from  methanol 
shortage  problem.  As  the  crossover  current  density  is  a  strong 
function  of  the  methanol  concentration,  the  crossover 
current  density  shown  in  Fig.  17  has  similar  distribution  trend  with 
Fig.  16. 


3.2.2.  Current  density  distribution 

As  shown  in  the  single  channel  model,  current  density  distri¬ 
bution  becomes  more  uniform  when  the  stoichiometry  increases 
(see  Fig.  18).  For  the  same  anode/cathode  stoichiometry  condition 
with  counter-flow  configuration,  the  highest  current  density  is 


Fig.  17.  Crossover  current  density  (A  nr2)  distribution  of  the  full-scale  cell  (28  cm2)  when  T„,u  =  313  K,  RH  =  10%,  CMeoH  =  4  M,  and  I  =  150  mA  cm“2:  (a)  la  =  2.0  and  lc  =  2.0  (case 
A),  (b)  la  =  1.4  and  =  2.0  (case  B),  (c)  la  =  2.0  and  £c  =  5.0  (case  C). 
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found  on  the  cathode  side  (low  half  region  in  case  A).  This  is 
because  ORR  is  more  active  than  MOR  and  the  anode  inlet  side 
suffers  from  methanol  crossover  which  degrades  the  cell 
performance. 


Note  cathode  inlet  port  (B  in  Fig.  18(a))  and  outlet  port  region 
(A  in  Fig.  18(a))  show  relatively  lower  current  density  than  other 
region  in  when  anode  stoichiometry  is  large  enough  (case  A  and 
C).  Firstly,  those  regions  are  corners  of  the  anode  channel  where 


a 


b 


c 


Fig.  18.  Current  density  (A  m~2)  distribution  in  the  membrane  of  the  full-scale  cell  (28  cm2)  when  Tan  =  313  K,  RH  =  10%,  CMboh  =  4  M,  and  /  =  150  mA  cnr2:  (a)  =  2.0  and 

=  2.0  (case  A),  (b)  =  1.4  and  =  2.0  (case  B),  (c)  =  2.0  and  ?c  =  5.0  (case  C). 
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fuel  concentration  is  greater  than  neighboring  region  due  to  the 
improved  mixing  effect  (see  Fig.  16).  Secondly,  those  regions  have 
the  widest  cathode  channel  area  (three  times  large  area  than 
a  single  channel  area)  since  two  cathode  channels  merges 
together  there,  which  means  oxygen  is  widely  spread  there. 
Therefore,  large  amount  of  methanol  crossovered  to  the  cathode 


side  reacts  very  actively  with  oxygen  there,  which  degrades  the 
cell  performance. 

3.2.3.  Temperature  distribution 

As  the  cathode  catalyst  generates  more  heat  that  the  anode 
catalyst,  temperature  distribution  follows  the  cathode  channel 
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shape  as  shown  in  Fig.  19.  Left  half  region  is  hotter  than  right  half 
region  since  methanol  crossover  concentrated  in  the  left  half  re¬ 
gion  generates  additional  heat.  Again,  hot  spots  are  found  in  region 
A  and  B  where  methanol  crossover  is  strong  and  cathode  area  is 
wide. 

3.2.4.  Effect  of  stoichiometry 

Comparing  case  A  and  case  B,  it  is  observed  that  increasing  fuel 
stoichiometry  improves  cell  voltage  while  worsening  fuel  efficiency 
(large  methanol  crossover).  Comparing  case  A  and  case  C,  it  is 
noticeable  that  increasing  air  stoichiometry  worsens  fuel  efficiency 
also  although  cell  voltage  improves.  This  is  because  improved  ox¬ 
ygen  reduction  reaction  in  the  cathode  catalyst  layer  due  to  large  air 
stoichiometry  induces  more  methanol  crossover  through  the 
membrane. 

4.  Conclusion 

Through  the  developed  non-isothermal  multi-dimensional 
DMFC  model  which  is  capable  of  simulating  saturation  jump  at 
the  interface  of  different  porous  media,  the  effectiveness  of  the 
both  anode  and  cathode  MPL  is  confirmed.  Non-isothermal 
behavior  including  thermal  diffusion  and  phase  change  is  exam¬ 
ined  and  it  is  found  that  heat  and  water  management  are  strongly 
coupled  with  each  other  by  the  material  property  (heat  conduc¬ 
tivity)  and  the  cell  operating  condition  (temperature  and  inlet  air 
humidity).  Finally,  local  distribution  of  current  and  methanol 
concentration  are  examined  with  a  full-scale  simulation  of  28  cm2 
DMFC  cell. 

List  of  symbols 


A  area  (m2) 

a  specific  reaction  area  (m2  m  3) 

ABP  ( abp )  anode  bipolar  plate 
ACL  ( acl )  anode  catalyst  layer 
ADL  ( adl )  anode  gas  diffusion  layer 
AML  (ami)  anode  micro-porous  layer 
cp  specific  heat  (J  l<g_1  K) 

c  concentration  (mol  m-3) 

CBP  (cbp) cathode  bipolar  plate 
CCL  ( ccl )  cathode  catalyst  layer 
CDL  ( cdl )  cathode  gas  diffusion  layer 
CML  ( cml )  cathode  micro-porous  layer 
d  Thickness  (m) 

D  diffusivity  (m2  s-1) 

F  Faraday  constant  (C  mor1) 

GDL  gas  diffusion  layer 

h  heat  transfer  coefficient  (W  m-2  K),  mass  transfer 

coefficient  (ms-1) 

I  current  (A) 

i  current  density  (Am-2) 

j  phase-diffusion  flux  (mol  m-2  s),  volumetric  current 

density  (A  m-3) 

J  Leverette  function 

k  thermal  conductivity  (W  m-1  K) 

k/i  Flenry’s  constant 

kr  relative  permeability 

Kc  rate  constant 

I<  permeability  (m2) 

m  mass  (kg) 

M  molar  weight 

MFM  (mem)  membrane 
MCO  methanol  crossover  rate 


MPL  (mpl)  micro-porous  layer 
N  molar  flux  (mol  m-2  s) 

p  pressure  (Pa) 

R  gas  constant  (J  K-1  mol) 

RH  relative  humidity  (%) 

s  liquid  saturation,  entropy 

S  source  term  (mol  m-3  s),  (W  m-3) 

T  temperature  (K) 

u  velocity  (ms-1) 

V  volume  (m3) 

Y  species  mass  fraction 

Greek 

a  transfer  coefficient,  net  water  transfer  coefficient 

5  Thickness  (m) 

e  porosity 

#  phase  potential  (V) 

p  overpotential  (V) 

k  ionic  conductivity  (S  m-1) 

A  relative  mobility,  water  content 

p  viscosity  (kg  m-1  s) 

v  kinetic  viscosity  (m2  s-1) 

6  angle  (°) 

p  density  (kg  m-3) 

a  Surface  tension  (N  m-1) 

?  stoichiometry 

r  shear  stress  (Pa) 

Q  ohmic 


Superscripts  and  subscripts 
A  anode 

b  breakthrough 

c  cathode,  contact 

capill  capillary 

d  drag 

diff  diffusion 

e  electrolyte  phase 

eff  effective 

g  gas 

hp  hydraulic  permeation 

i  species 

k  species 

1  liquid 

m  mass 

mem  membrane,  methanol 

0  reference 

ref  reference 

s  solid  phase 

sat  saturated 

tot  total 

T  temperature 

xover  crossover 
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